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Abstract 

The density matrix renormalization group method is used to investigate the spin-Peierls tran- 



polates between a gapped, dispersionless (Einstein) limit to a gapless, dispersive (Debye) limit. A 



sition for Heisenberg spins coupled to quantized phonons. We use a phonon spectrum that inter- 
U 
CO 

variety of theoretical probes are used to determine the quantum phase transition, including en- 

8 . 

ergy gap crossing, a finite size scaling analysis, bond order auto-correlation functions, and bipartite 
quantum entanglement. All these probes indicate that in the antiadiabatic phonon limit a quantum 

o' 

O . phase transition of the Berezinskii-Kosterlitz-Thouless type is observed at a non-zero spin-phonon 

P^ | coupling, g c . An extrapolation from the Einstein limit to the Debye limit is accompanied by an 

> ■ 

f^*) ' increase in g c for a fixed optical (q = n) phonon gap. We therefore conclude that the dimerized 

ground state is more unstable with respect to Debye phonons, with the introduction of phonon 
i^ , dispersion renormalizing the effective spin-lattice coupling for the Peierls- active mode. We also 

o 

show that the staggered spin-spin and phonon displacement order parameters are unreliable means 
of determining the transition. 
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I. INTRODUCTION 

Since the discovery of high-temperature superconductivity in doped antiferromagnets 
there has been a marked increase in interest - both theoretical and experimental - in low- 
dimensional quantum magnetism. However, the effect of the interaction of quantum spins 
with further degrees of freedom such as disorder, phonons, and holes produced by doping 
remains relatively poorly understood. 

The instability of the spin-| antiferromagnetic Heisenberg chain to a static uniform distor- 
tion gives rise to the so-called spin-Peierls (SP) transition. This occurs because the explicit 
dimerization opens a gap, A, in the spin-excitation spectrum, lowering the total magnetic 
energy by an amount that offsets the accompanying increase in lattice energy. 

The SP instability is itself the antiferromagnetic analogue of the Peierls transition^: a half- 
filled one- dimensional metallic phase is unstable with respect to a commensurate periodic 
lattice distortion of wave vector q = 2k p. Indeed, SP models represent the large on-site 
coupling limit of the corresponding half-filled Hubbard-Peierls Hamiltonian: for infinite 
inter-electron repulsion, charge degrees of freedom are effectively quenched resulting in the 
loss of electron itinerancy for a half-filled band. For linear chains under open boundary 
conditions the Jordan- Wigner transformation maps the Heisenberg-SP model onto a spinless 
fermion- Peierls model with nearest-neighbour repulsion, highlighting the decoupling of spin 
and charge degrees of freedom and explaining the SP nomenclature. 

The SP instability is well understood in the static-lattice limit for which the frequency, 
uv, of the Peierls- active mode is taken to be much smaller than the antiferromagnetic ex- 
change integral, J. In this adiabatic phonon limit the ground state (GS) is known to have 
a broken-symmetry staggered dimerization for arbitrary electron-phonon (e-ph) coupling. 
Experimentally, such behavior was first observed in the 1970s for the organic compounds of 
the TTF and TCNQ series^. For many quasi-one-dimensional materials, however, the zero- 
point fluctuations of the (quantized) phonon field are comparable to the amplitude of the 
Peierls distortion^-. In CuGeOs 6 , for example, Cu 2+ ions form well separated spin-| chains 
with an exchange interaction that couples to high-frequency optical phonons w T ~ ®{J)- 
CuGe0 3 has since become a paradigm of inorganic antiadiabatic SP behaviour, stimulating 
several numerical studies of dynamical phonon models^-. 

Using both the density-matrix renormalization group (DMRG) method^ and renormal- 



ization group (RG) methods within a bosonization scheme^, it has been demonstrated that 
quantum fluctuations destroy the Peierls state for small, non-zero couplings in both the spin- 
less and spin-| Holstein models at half-filling. Analogous results for the XY-SP model with 
gapped, dispersionless (Einstein) phonons were also obtained by Caron and Moukouri^, us- 
ing finite-size scaling analysis of the spin gap to demonstrate a power-law relating the critical 
coupling and the Peierls-active phonon frequency: g* Y ~ uj®' 7 . For models with sufficiently 
large Einstein frequency, gapped phonon degrees of freedom can be integrated away to gen- 
erate a low-energy effective-fermion Hamiltonian characterized by instantaneous, non-local 
interactions^. For spinless models, RG equations indicate that unless the non-local contri- 
bution to the umklapp term has both the right sign and a bare (initial) value larger than a 
certain threshold, the umklapp processes are irrelevant and the quantum system is gapless 13 . 
Conversely, if the threshold condition is satisfied, the umklapp processes and vertex function 
grow to infinity, signalling the onset of gapped excitations and a dimerized lattice. 

For the Su-Schrieffer-Heeger (SSH) model, Fradkin and Hirsch undertook an extensive 
study of spin-i (n = 2) and spinless [n = 1) fermions using world-line Monte Carlo 
simulations^. In the antiadiabatic limit (i.e. vanishing ionic mass M), they mapped the 
system onto an n-component Gross-Neveu model, known to exhibit long-ranged dimeriza- 
tion for arbitrary coupling for n > 2 (although not for n — 1). For M > an RG analysis 
shows the low-energy behavior of the n = 2 model to be governed by the zero-mass limit of 
the theory, indicating that the spinful model presents a dimerized GS for arbitrarily weak 
e-ph couplings. The spinless model, on the other hand, has a disordered phase for small 
coupling if M is finite, with an ordered phase realized for bare coupling in excess of a certain 
threshold. As M — > oo the size of the disordered region shrinks to zero, reconnecting with the 
adiabatic result of Peierls and Frolich 1 . Later work by Zimanyi et alM on one-dimensional 
models with both electron-electron (e-e) and e-ph interactions showed they were found to 
develop a spin gap if the combined backscattering amplitude gj = giius) +gi{uj) < 0, where 
gi{uj) is the contribution from electron-electron (e-e) interactions and gi(oo) < is the e- 
ph contribution in the notation of~. Hence, for the pure spinful SSH model, g\ = and 
gj < for any nonzero e-ph coupling, implying a Peierls GS for arbitrary e-ph coupling, in 
agreement with the earlier MC results^. 

In this article we study the influence of gapless, dispersive antiadiabatic phonons on the 
GS of the Heisenberg-Peierls chain. That this model is yet to receive the same level of 



attention as its gapful, dispersionless counterpart is due in part to the presence of hydro- 
dynamic modes, resulting in logarithmically increasing vibrational amplitudes with chain 
length. To this end, the authors of-^ ancr^ assumed acoustic phonons to decouple from the 
low-energy spin states involved in the SP instability, motivating the retention of only the 
optical phonons close to q = tt. In this regard, optical phonons have been expected to be 
equivalent to fully quantum mechanical SSH phonons. For pure Einstein phonons, Wellein, 
Fehske, and Kampf-, however, found that the singlet-triplet excitation is strongly renormal- 
ized when phonons of all wavenumber are taken into account, the restriction to solely the 
q = tc modes leading to a substantial overestimation of the spin gap. Physically, this implies 
that the spin-triplet excitation is accompanied by a local distortion of the lattice, neces- 
sitating a multiphonon mode treatment of the lattice degrees of freedom. We anticipate, 
then, that truncating the Debye-phonon spectrum to leave only those modes which couple 
directly to the SP phase may be not be physically reasonable. 

In this work we use the DMRG technique to numerically solve the Heisenberg-Peierls 
model with a generalized gapped, dispersive phonon spectrum. The phonon spectrum inter- 
polates between a gapped, dispersionless (Einstein) limit and a gapless, dispersive (Debye) 
limit. We proceed by considering a system of Heisenberg spins dressed with pure Einstein 
phonons for which we observe a Berezinskii-Kosterlitz-Thouless (BKT) quantum phase tran- 
sition at a non-zero spin-lattice coupling. Progressively increasing the Debye character of the 
phonon dispersion (at given phonon adiabaticity) results in an increase in the critical value 
of the spin-lattice coupling, with the transition remaining in the BKT universality class (see 
Section IIIICI) . These findings are corroborated by an array of independent verifications: 
energy-gap crossings in the spin-excitation spectra (see Section [HI Aj) . finite-size scaling of 
the spin-gap (see Section [IIIBI) . bond order auto-correlation functions (see Section FlIIDj) . 
and quantum bipartite entanglement (see Section UlIEj) . 



We note that earlier DMRG investigations of the Heisenberg-SP Hamiltonian with Debye 
phonons indicated a dimerized GS for arbitrary coupling^. This conclusion was based on 
the behavior of the staggered phonon order parameter, m p , (defined in Section IIIIDI) . In 
this paper we show that m p is an unreliable signature of the transition. 

In the next Section we describe the model, before discussing our results in Section IIHI 



II. THE MODEL 

The Heisenberg spin-Peierls Hamiltonian is defined by, 

H = H s . p + H p . (1) 

H s _ p describes the spin degrees of freedom and the spin-phonon coupling, 

H s _ p = ^2\J + a(ui +1 - ui)]Si ■ Si +1 , (2) 

i 

where S; is the Pauli spin operator, ui is the displacement of the Zth ion from equilibrium, 

and a is the spin-phonon coupling parameter. 

H p describes the lattice degrees of freedom. In the Einstein model the ions are decoupled, 

H E = yEL + \ K y u i (3 ) 

p ^ 2M 2 ^ 

i i 

In the Debye model, however, the ions are coupled to nearest neighbors, 

i i 

For the Einstein phonons it is convenient to introduce phonon creation, bj, and annihila- 
tion operators, 6/, for the Zth site via, 

"-(s^) 1/2(6!+M (5) 

and 

.{ Mhcu x \ 1/2 ^ 
Pl = l \~~2 — ) ( b i~ bl )> ( 6 ) 

where 

u x = uj e = ^/K/M = u b . (7) 

Making these substitutions in Eq. (j2]) and Eq. (J3J) gives, 



H s . p = J 22 

i 

and 



H-^(^) 1/2 (B,-B J+1 ) 



< = ^E(W£) 



S L ■ S m (8) 



(9) 



where Bi = \{b\ +bi) is the dimensionless phonon displacement and, 



9e = ol 



Mu 2 E J 



1/2 



a 



2 , 



\ 1/2 



KJ 






(10) 



is the dimensionless spin-phonon coupling parameter. 

For the Debye phonons we introduce phonon creation and annihilation operators defined 
by Eq. ([5]) and Eq. (jSJ) where 



UJ X 



uj d = ^2K/M = V2oo b . 



(11) 



Making these substitutions in Eq. (J2J) and Eq. (J3J) gives, 



H s _ p = J 22 



hup X 



i+9D[-f) [B l -B l + X ] 



S; ■ S/ + i 



and 



where, 



H° = hw D J2 Uk + ^J-hw D J2 B\ +1 B U 

1 \ 1/2 



1/2 



gD = a \M^j) =a \KJ 






H^ may be diagonalized by a Bogoluibov transformation to yield 



H° = hJ2u D (q)P% 



(12) 



(13) 



(14) 



(15) 



where ojd(q) is the dispersive, gap less phonon spectrum, 

LO D (q) = V2uJd sin (-) , (16) 

for phonons of wavevector q. 

We now introduce a generalized spin-phonon model with a dispersive, gapped phonon 
spectrum, via 



H s _ p = J 22 



I + 9 [ ^j) 1/2 (B l - B l+1 ) 



Si ■ Si + i 



and 



H p = h(u E + u D ) J2 \ h \ h i + o ) " fk0D J2 B w Bl 
Again, Eq. (|T5|) may be diagonalized to give, 

H p = hy^u(q)/3g/3 q + constant, 



(17) 



(18) 



(19) 



where, 

u(q) = (u E + w fl ) ( 1 - ( U ° ) cos q] , (20) 

is the generalized phonon dispersion, as shown in Fig. [TJ 
The q = phonon gap frequency is, 

u(q = 0) = tu = (ue(ue + wb)) 1/2 (21) 

and the q = n optical phonon frequency is, 

u(q = 7r) = ^ = ((2w B + u D )(u E + w D )) 1/2 . (22) 

We now define the dispersion parameter 7 as, 

7 = w /W (23) 

7 is a mathematical device that interpolates the generalized model between the Einstein 
(7 = 1) and Debye (7 = 0) limits for a fixed value of the q = i\ phonon frequency, u n . 
The dimensionless spin-phonon coupling, g, as well as u n /J and 7 are the independent 
parameters in this model, ue and u D , on the other hand, are determined by Eq. (I2T|) . ff22|) . 
and ((23]). 

The generalized model can be mapped onto the Einstein and Debye models by the ob- 
servation that in the Einstein limit, 



while in the Debye limit, 



co n = cue = co b = y/K/M; (24) 

9 = 9e, 

uj w = v2cud = 2ub, (25) 

9 = W2 1/4 - 

The introduction of a generalized phonon Hamiltonian avoids the problems associated 
with hydrodynamic modes and places a criterion on the reliability of the gap-crossing char- 
acterization of the critical coupling (as described in Section IIII Al) . Starting from the 



Heisenberg-SP Hamiltonian in the Einstein limit (7 = 1), the effect of dispersive lattice 
fluctuations can be investigated via a variation of 7. The Debye limit is then found via an 
extrapolation of 7 — > 0. 



.3 




FIG. 1: Generalized phonon dispersion, w(g), denned in Eq. (f20|) . (1 — j)co n is the phonon 'band 
width' (which vanishes in the Einstein-limit), while 70% is the phonon 'mass-gap' (which vanishes 
in the Debye-limit). The dispersion parameter, 7, and the optical phonon frequency, uj n , are model 
parameters. 

The model is solved using the density matrix renormalization group (DMRG) method 18 
with periodic boundary conditions throughout. Our implementation of the DMRG method, 
including a description of the adaptation of the spin-phonon basis and convergence, is given 
in the Appendix. 

III. RESULTS AND DISCUSSION 
A. Gap-crossing 



For the Einstein model with a non-vanishing value of Ue the critical spin-phonon cou- 
pling, g c , may be determined using the gap-crossing method of Okamoto and Nomura^, 
as shown in Fig. [2] for an 80-site chain. If the iV-site system has quasi-long-range Neel 
order for < g < g c (N), the lowest excitation is a triplet state, i.e. A st < A ss and 



lim 



N-yoo A st 



lim 



iV— >oo ^ss 



0, where A st and A ss are the triplet and singlet gaps, re- 
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FIG. 2: Gap-crossing construction for the 7 = 1 (Einstein) Heisenberg-SP model for N 
Inset: The infinite-chain critical coupling g£° is determined by extrapolation. 
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spectively. Conversely, for g > g c (N), the system is dimerized with a doubly-degenerate 
singlet GS in the asymptotic limit (corresponding to the translationally equivalent 'A' and 
'B' phases), while the lowest energy triplet excitation is gapped. However, for finite systems 
the two equivalent dimerization phases mix via quantum tunneling, and now A ss < A st , 
with limAr^oo A ss = and liniAr^oo A st > 0. The gap-crossing condition A st = A ss therefore 
defines the finite-lattice crossover coupling g c (N). 

For the Debye model, however, the gap-crossing method fails because of the q — > 
phonons that form a gapless vibronic progression with the groundstate. The hybrid spectrum 
(shown in Fig. [1]) allows us to extrapolate from the pure Einstein limit to the Debye limit, as 
the lowest vibronic excitation is necessarily ju^. Provided that A ss < u(q = 0) = ju^, the 
gap crossover method unambiguously determines the nature of the GS. We can confidently 
investigate Eq. (0Q) for (0.1 < 7 < 1) with u n /J E [1, 10], thereby determining g c (iV, 7). A 
polynomial extrapolation of 1/N — > generates the bulk-limit critical coupling <j£° for a 




FIG. 3: Phase diagram in the g^-uo^ plane for the infinite Heisenberg-SP chain for 7 = 1 (diamonds) 
and 7 = 0.5 (triangles); extrapolation to 7 = generates the Debye-limit (open circles). Inset: 
variation of g^ with 7 for the antiadiabatic limit, w f / J = 10. 

given 7 (as illustrated in Fig. [2]). A subsequent polynomial extrapolation determines the 
7 = (Debye) limit. A phase diagram for the Heisenberg-SP chain found in this way is 
shown in Fig. |3j Notice that for a fixed u n the critical coupling is larger for the Debye 
model than for the Einstein model, showing that the quantum fluctuations from the q < it 
phonons (as well as the q = n phonon) destablize the Peierls state. 

Following Caron and Moukouri 11, we tentatively propose a general power-law for the 
Heisenberg-SP model, relating the bulk-limit critical coupling to u n for a given 7, 



5c »k,7) = /3(7M^. 



(26) 



The infinite-chain values of /3 and 77, and g™ for u^/ J = 10 are given in Table HI We find a 
non-zero critical coupling for all phonon regimes 7, with the absolute value of g%° increasing 
as 7 — > 0, as shown in the inset of Fig. |3j 
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0.511 


0.437 


1.397 
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0.452 


0.392 


1.103 


1 (Einstein) 


0.350 


0.337 


0.761 



TABLE I: Gap-crossing determined bulk-limit values of /3(7) and 77(7) (denned by Eq. (|26p ). and 
g£° for uj n /J = 10. The Debye limit (7 = 0) is obtained by extrapolation of 7 — > 0. 

B. Finite-size scaling 



In order to ascertain the analytic behavior of the spin gap from the numerical data 
it is necessary to account for finite-size effects. We assume that the (singlet-triplet) gap 
Ajv = A st for a finite system of N sites obeys the finite-size scaling hypothesis^^i 

1 



A N = ^F(NA C 



(27) 



with Ago the spin-gap in the bulk limit. Recalling that g%° = lim A r_ Sh00 f7 c (A^), it follows that 
Aoo(Sc°) = and so curves of NA^ versus g are expected to coincide at the critical point 
where the bulk-limit spin-gap vanishes, as confirmed in Fig. H] and Fig. [51 

The finite-size scaling method is more robust than the gap-crossing approach, being 
applicable to the SP Hamiltonian for all values of 7. On the other hand, its use as a 
quantitative method is limited by the accuracy with which plots may be fitted to Eq. (1271) . 
In practice, plots of NA st (N) versus g become progressively more kinked about the critical 
point as 7 — > 0. Nevertheless, we find F to be well approximated by a rational function and 
the resulting g£°(7) to be in accord with the predictions of the gap-crossover method. 



C. Berezinskii-Kosterlitz-Thouless transition 



For a BKT transition the spin-gap A = hin^^oo A st is expected to exhibit an essential 
singularity at <j£° with plots of A st versus g for iV — > 00 found to be well fitted by the Baxter 
forrn^ (as shown in Fig. [6]), 

A~af(g)exp(-b[f(g)} 2 ) (28) 

where^, 

f(g) = (9-9?)- 1/2 - (29) 
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FIG. 4: NA st (N) versus the spin-phonon coupling, g, for the 7 = 0.5 Heisenberg-SP model for 
N = 16 (diamonds), 40 (squares), 80 (triangles), and 160 (open circles) for uj n /J = 10. The curves 
converge at g£° (the value shown is obtained via gap-crossing). 

Extrapolating A st (N) for 1/iV — > generates A for a given 7 and it is possible, in 
principle, to distinguish dimerized from spin-fluid GSs by examining the scaling behavior 
of Agt(iV), which tends to zero in the bulk- limit for the spin fluid and to a non-zero A for 
the gapped phase. However, not only must three parameters (a, b, and g%°) be obtained 
from a non-linear fit (shown in Table [XT]) , but there is considerable difficulty in determining 
A accurately near the critical point: the spin-gap is extremely small even for values of g 
substantially higher than <j£° due to the essential singularity in Eq. (1281) . Determining such 
small gaps from finite-size scaling is highly problematic with very large lattices required to 
observe the crossover from the initial algebraic scaling (in the critical regime) to exponential 
scaling (for gapped systems). Hence the gap-crossover method is expected to be substantially 
more accurate than a fitting procedure for the determination of the critical coupling, the 
latter tending to overestimate <j£° (see ref^), as confirmed by a comparison of Tables [land 
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FIG. 5: NA st (N) versus the spin-phonon coupling, g, for the 7 = (Debye) Heisenberg-SP model 
for N = 16 (diamonds), 40 (squares), 80 (triangles), and 160 (open circles) for co n /J = 10. The 
curves converge at g c (the value shown is obtained via gap-crossing). 
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4.110 
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14.206 


3.042 


0.731 



TABLE II: Baxter-equation parameters obtained by fits to Eq. (128p for u} n /J = 10. 



D. Correlation functions and order parameters 



The q = 7T structure factor, S(q), of the bond-order auto-correlation function can be used 
to determine the phase transition. S(q) is defined by, 



S(q) = \exp(igm)(C(m)) 



(30) 



where 



^M = 7V E(°< " (°l))(°l+™ ~ <Ol+m». 



(31) 
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FIG. 6: Bulk-limit singlet-triplet gap, A, as a function of the spin-phonon coupling, g, with 7=1 
(Einstein) [diamonds] 7 = 0.5 [open squares] and 7 = (Debye) [open circles] for lj w / J = 10. Plots 
are fitted to the BKT form (Eq. ([28]) ). 



The 'bond order' operator, Oi, is given by 



O, 



S*Sf +1 spin-spin; 

(qi — qi+x) phonon displacement. 



The transition to a dimerized state is marked by the development of a staggered kinetic 
energy modulation and quasi long-range-order in the bond order. It is signalled by a di- 
vergent peak in S(q = tt) at the critical coupling in the asymptotic limit, as shown in Fig. 
[3 (We note, however, that the structure factor associated with the phonon displacement 
bond order auto-correlator, 0\ = (qi — qi+i), fails to resolve the transition, instead increasing 
monotonically with g for all 7.) 

Followingi£>22 we also consider staggered spin-spin, m s , and phonon displacement, m p , 
order parameters, 

^^^E^^^m) ( 32 ) 
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FIG. 7: Structure factor for the spin-spin bond order auto-correlation function at q = ir for the 
7 = (Debye) model with uj w /J = 10 for N=20 (squares), iV=40 (triangles), and iV=80 (open 
circles). 



and 



m T 



^£(-l)'<5m-^>. 



(33) 



For linear chains under OBC, the end sites break the energetic degeneracy between the 
otherwise equivalent \A) and \B) states, which are related by a translation of one repeat 
unit. Physically, however, PBC are strongly preferable to OBC as boundary effects are 
eliminated and finite-size extrapolations can be performed for much smaller N. In addition, 
greater accuracy can be obtained by investigating cyclic chains, although it is necessary 
to explicitly break the degeneracy between the 'A' and 'B' phases through inclusion of a 
symmetry-breaking term H' in Eq. (pQ), 



ff' = p^(-l)'(S r S w ), 



(34) 



and extrapolating p — y 0. 

Our results indicate that both m s and m p scale to zero as g — > oo, suggesting g£° = + for 
all 7, accounting for the earlier findings of 16 . These predictions are incomplete deviation to 
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the other positive signatures of a phase transition for g > 0. We attribute this discrepancy to 
the action of the perturbation H' (Eq. (jM|) ) on the fixed-point behavior of the Heisenberg-SP 
Hamiltonian for an insufficiently small perturbation, and thus conclude that the staggered 
order parameters must be treated with caution when determining the phase transition. 

E. Quantum bipartite entanglement 

It has recently been conjectured that quantum entanglement plays an important role in 
the quantum phase transitions (QPT) of interacting quantum lattices. At the critical point- 
as in a conventional thermal phase transition — long-range correlations pervade the system. 
However, because the system is at T = (and assuming no ground-state degeneracy) the 
GS is necessarily a pure state. It follows, then, that the onset of (long-range) correlations - 
being the principal experimental signature of a QPT - is due to entanglement in the GS on 
all length scales. 

For an iV-site lattice, bipartite entanglement is quantified through the von Neumann 



entropy 24 , 



Sl = -Trsp<j(L) log 2 p s (L) = - ^ "* log 2 v a , (35) 



where ps(L) is the reduced-density matrix of an L-site block (typically coupled to an L-site 
environment S such that 2L = N) and the u a are the eigenvalues of ps{L). It is clear 
from Eq. fl35l) that a slow decay of the reduced density-matrix eigenvalues corresponds to a 
large block entropy. Provided the entanglement is not too great and the v a decay rapidly, a 
matrix-product state is then a good approximation to the GS 2 ^. We note here the utility of 
the DMRG prescription in determining Sl = S^/t^- 

Wu et al. 27 argued, quite generally, that QPTs are signalled by a discontinuity in some 
entanglement measure of the infinite quantum system. For finite one-dimensional gapped 
systems the noncritical entanglement is characterized by the saturation of the von Neumann 
entropy with increasing L: the entropy of entanglement either (i) vanishes for all L or 
(ii) grows monotonically with L until it reaches a saturation value for some block length 
Lq 28 . Noncritical entanglement in the GS corresponds, thus, to a weak, semi-local 29 form of 
entanglement driven by the appearance of a length scale L due to, e.g, a mass gap in the 
Hamiltonian. For any L, the reduced-density matrix ps(L) is effectively supported on just 
a small, bounded subspace of the L-spin Hilbert space. Critical models, on the other hand, 
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FIG. 8: Von Neumann entropy, Sl, for the 7 = (Debye) model with u> n /J = 10 for lattice sizes 
N = 20 (squares), 40 (triangles), and 80 (open circles); L = N/2. 

are expected to exhibit logarithmic divergence in Sl at large L: Sl = fclog 2 L + constant. 

For a given total system size N and phonon dispersion 7, the block entropy is found to be 
maximal for a non-zero spin-phonon coupling g c (N), close to the corresponding gap-crossing 
and bond order structure factor values (as shown in Table IIIII) . As shown in Fig. [SJ in 
the critical regime, g < g c (N), the block entropy is indeed found to scale logarithmically 
with system-block length, while in the gapped phase, g > g c (N), it is characterized by 
the emergence of a saturation length scale L$ that varies with 7. These findings are in 
agreement with 28 and consistent with the observation that the transition belongs to the 
BKT universality class^. 



F. Phase diagram 

To conclude this section we discuss the phase diagram of the Heisenberg spin-Peierls 
model. Fig. [3] shows the phase diagram as a function of the model parameters g and the 
q = it phonon gap, u> n , as defined in Eq. (TT7|) and Eq. (TT9|) . Evidently, for a fixed value of w T 
the spin-Peierls state is less stable to dispersive, gapless quantum lattice fluctuations than 
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N 


gap 


flf 


<7c vN 


20 


1.441 


1.441 


1.443 


40 


1.419 


1.420 


1.422 


80 


1.400 


1.402 


1.403 



TABLE III: Consistency of the various probes of the transition: critical spin-phonon couplings 
determined by gap-crossing (gap), q = it structure factor of the bond order auto-correlation function 
(SF), and von Neumann entropy (vN) for N = 20, 40, and 80 sites. 7 = (Debye) and oj^/J = 10. 

to gapped, non- dispersive fluctuations, implying that the q < tt phonons also destablize the 
Peierls state. 

It is also instructive, however, to plot the phase diagram as a function of the physical 
parameters a and ujb = y/K/M, as defined in Eq. (J2J), Eq. (j3]) and Eq. (@J. The mapping 
between model and physical parameters is achieved via Eq. (TTUj) . Eq. (fl4"l) . Eq. ( 124"|) . and 
Eq. fl25|) (and setting K — 1). Since u> n = ujb for the Einstein model, whereas u n = 2ub 
for the Debye model, the Debye model is further into the antiadiabatic regime for a fixed 
value of ujb- We also note that for a given model electron-phonon coupling parameter, g, 
the physical electron-phonon coupling parameter, a, is larger in the Debye model than the 
Einstein model (see Eq. ( ITU]) and Eq. (1T4l)). Consequently, we expect the dimerized phase 
to be less robust to quantum fluctuations in the Debye model for fixed values of Ub and a, 
as confirmed by Fig. HP 1 . 

IV. CONCLUSIONS 

The coupling of spin and lattice degrees of freedom with reduced dimensionality results 
in the instability of a one-dimensional Luttinger liquid towards lattice dimerization and the 
opening of a gap at the Fermi surface. Coupling to the lattice gives rise to a BKT transition 
from a spin liquid with gapless spinon excitations to a dimerized phase characterized by an 
excitation gap. 

For the quantum Heisenberg chain in the antiadiabatic limit (J/u n << 1), the spin- 
fluid phase becomes unstable with respect to lattice dimerization above a non-zero spin- 
phonon-coupling threshold for all phonon gaps, joo^. This observation holds for u n /J ~ 1. 
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FIG. 9: Phase diagram in the a£° -a;& plane for the infinite Heisenberg-SP chain for 7=1 (Einstein), 
diamonds; and 7 = (Debye), squares. 



Increasing the contribution of dispersive phonons to H p gives rise to an increase in the critical 
coupling, supporting the intuition that gapless phonons more readily penetrate the GS (with 
the q < 7T phonon modes renormalizing the dispersion at the Peierls-active modes). This 
observation has been corroborated by an array of independent verifications. The behavior of 
the Debye model is qualitatively different from the Einstein model, with different exponents 
and prefactors for the critical coupling versus phonon frequency, Eq. ( 126|) . and the Baxter 
expression for the spin- gap, Eq. ( l28l) . 

We note that staggered order parameters are an unreliable means of determining the 
phase transition^, because of the use of a symmetry breaking perturbation for PBCs that 
changes the FP of the Hamiltonian for insufficiently small perturbations. 

Placing these findings in the context of experiment, estimates of the model parameters for 
a number of spin-Peierls compounds are listed in Table HVl (reproduced from 9 ). Clearly, the 
static approximation is not applicable to CuGeOs. In addition, it is also questionable as to 
whether the listed organic spin-Peierls compounds are themselves truly static lattice mate- 
rials, thereby justifying a dynamical phonon treatment. The most physically relevant region 
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of the phonon spectrum, however, appears to be one of intermediate frequency, dividing the 
adiabatic and antiadibatic limits. Nevertheless, even though cj^/J = 3 for CuGeC>3, refer- 
ring to the phase diagram of Fig. [3] we note the applicability of the generalized power-law for 
small u) n /J, and hence the occurrence of a finite critical coupling in the regime applicable 
to CuGe0 3 . 



Material 


J 


w,r 


A 


CuGe0 3 

TTFCuS 4 C4(C 3 F)4 
(MEM)(TCNQ) 2 


100 

70 
50 


300 
10* 
100 


20 
20 
60 



TABLE IV: Estimates of the antiferromagnetic exchange, J, dimerization phonon frequency, u n , 
and spin gap, A, for various SP materials. All units are in Kelvins. * Value deduced by comparison 
of experimental values of J, the transition temperature, and A to mean-field theoretical expressions 
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Examination of intermediate phonon frequencies and their extension to spinful fermion 
models with Coulomb repulsion^ is straightforward and currently in progress. 

Appendix A: DMRG and in situ optimization 

We solve Eq. ( !T7|) and Eq. ( !T8|) using the real-space density matrix renormalization group 
(DMRG) method 18 , with ten oscillator levels per site, typically ~ 200 block states and ca. 
10 6 superblock states. Finite lattice sweeps are performed at target chain lengths under 
PBC. The convergence indicators are shown in Tables IVllIXt with additional convergence 
tables in ref-^ for the same model. 

The DMRG algorithm typically proceeds through the augmentation of an (L — l)-site 
system block S by a single site i a . The augmented system block 5" (L sites) is coupled 
to an augmented L-site environment block E', formed analogously to S'. The system and 
environment comprise the superblock (2L = N sites), whose state vector, |\1/), is readily 
obtained by a suitable diagonalization routine. By tracing over the degrees of freedom in 
E', the reduced-density matrix of S' (ps>{L) = Tr E i\^/)(^/\) is obtained; a pre-determined 
proportion of the largest-eigenvalue eigenstates of ps>{L) is retained, forming the system- 
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block basis for the next iteration. The DMRG prescription results in (9(iV)-growth of the 
superblock Hilbert space. 

In addition (and prior) to the effective truncation and rotation of the system-block basis, 
we employ a single-site optimization at each DMRG step. All superblock degrees of freedom, 
save those belonging to i a , are traced over and the resulting single-site reduced-density matrix 
p a is diagonalized, generating an optimal single-site basis in the correct physical environment 
with which to augment ,§LiSj22i22jM. This local Hilbert space adaptation generates single- 
site bases, the principal utility of which is the solution of many-body problems with large 
numbers of degrees of freedom. A controlled truncation of a large Hilbert space therefore 
allows a small (for our purposes six-dimensional) optimal basis to be used without significant 
loss of accuracy, making it ideally suited to spin-phonon problems, where the number of 
phonons is not conserved and the phonon Hilbert space is, in principle, infinite. 



m 


E g /J 


ni 


cr n 


2 


-18.746372 


0.00291 


0.0468 


5 


-18.747323 


0.00221 


0.0470 


8 


-18.747323 


0.00221 


0.0470 


10 


-18.747367 


0.00221 


0.0470 



TABLE V: GS energy, E g /J, average phonon occupation number, m = (b]bi), and standard 
deviation, a n , for a 40-site chain with 7 = 1 (Einstein), u w /J = 10, and m oscillator levels 
per site. 



m 


E g /J 


ni 


<Jn 


2 


-27.198068 


0.0331 


0.179 


5 


-31.795996 


0.1477 


0.437 


8 


-31.991408 


0.1717 


0.481 


10 


-31.991541 


0.1719 


0.484 



TABLE VI: GS energy, E g /J, average phonon occupation number, ni, and standard deviation, a n , 
for a 40-site chain with 7 = (Debye), uj n /J = 10, and m oscillator levels per site. 
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m 


E g /J 


n\ 


&n 


2 


-18.882263 


0.0318 


0.175 


5 


-19.358486 


0.1368 


0.420 


8 


-19.371279 


0.1520 


0.451 


10 


-19.371877 


0.1522 


0.452 



TABLE VII: GS energy, E g /J, average phonon occupation number, ni, and standard deviation, 
a n , for a 40-site chain with 7 = (Debye), cj^/J = 1, and m oscillator levels per site. 



e 


E g /J 


M 


SBHSS 


10- 10 


-24.266060 


176 


14572 


10- 11 


-24.350806 


278 


32328 


10- 12 


-24.377641 


300 


60054 


10- 15 


-24.377646 


350 


172654 



TABLE VIII: GS energy, E g /J, of Heisenberg-SP model as a function of the density-matrix eigen- 
value product cutoff, e, number of system block states, M, and the superblock Hilbert space size, 
(SBHSS) for a 40-site chain with 10 oscillator levels per site and 7 = 0.5. 
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